Determination of step— edge barriers to interlayer transport from surface morphology 
during the initial stages of homoepitaxial growth 
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We use analytic formulae obtained from a simple model of crystal growth by molecular-beam 
epitaxy to determine step-edge barriers to interlayer transport. The method is based on infor- 
mation about the surface morphology at the onset of nucleation on top of first-layer islands in 
the submonolayer coverage regime of homoepitaxial growth. The formulae are tested using kinetic 
Monte Carlo simulations of a solid-on-solid model and applied to estimate step-edge barriers from 
scanning-tunneling microscopy data on initial stages of Fe(OOl), Pt(lll), and Ag(lll) homoepitaxy. 
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Nearly thirty years ago, it was found in field-ion- 
microscopy experiments that interlayer transport on 
some metal surfaces is suppressed by additional acti- 
vation barriers to hopping over step edges, 1 and such 
barriers were recently shown to be present on a semi- 
conductor surface as well. 2 The consequences of this 
Ehrlich-Schwoebel (ES) barrier for growth on vicinal 3,4 
and singular 5 surfaces have been theoretically investi- 
gated. In particular, study of the effects of the ES bar- 
rier for growth on a singular surface led to new insights 
into homoepitaxial growth modes 6 as well as into ki- 
netic roughening of growing surfaces. 2 The ES barrier 
has emerged as a material parameter of an importance 
comparable to the surface diffusion barrier. It is, how- 
ever, very difficult to determine this quantity experimen- 
tally. 

Very recently, Meyer et al. 7 proposed (amongst other 
suggestions) a method of determining the ES barrier 
based on the surface morphology at the onset of nucle- 
ation on top of first-layer islands as seen by scanning 
tunneling microscopy (STM). In this paper, we provide 
a more consistent analytical treatment of this problem 
and derive a formula different from the one proposed 
by Meyer et al.. This formula and its modifications are 
tested using kinetic Monte Carlo (KMC) simulations and 
applied to STM data for Fe/Fe(001), Ag/Ag(lll), and 
Pt/Pt(lll) homoepitaxy. 

We consider a model of the initial stages of homoepi- 
taxy similar to the one proposed by Lewis and Anderson 8 
with circular-shaped islands regularly distributed over a 
perfect singular surface (Fig. 1). Growth is initiated by 
a flux of atoms incoming to the surface. The adatom 
density p on the surface increases until it reaches a crit- 
ical value p c at a time t = when islands of radius 7*0 = 1 
separated by a distance 2L are formed. (Notice that all 
lengths are given in units of the lattice constant and are 
accordingly dimensionless.) We assume that the inter- 
island free-adatom density is then well approximated by 




FIG. f. Model geometry 



its steady-state form with weakly time-dependent coef- 
ficients (quasi-static approximation 9 ). We will consider 
relaxing the assumption of an instantaneous steady state 
later. Each island has an effective catchment area of ra- 
dius L with free-adatom flux equal to zero at a distance 
r = L. Based on these assumptions, the adatom density 
on the substrate immediately following nucleation is the 
solution of the diffusion equation 

DV 2 p + F = (1) 

with boundary conditions (cf. Ref. 9, Section 7.4) 
dp/dr \ r =R(t) =p(R(t)) (and thus also dp/dr | r=ro =_R(o) = 
p(ro)) and dp/dr \ r= L = 0, where R(t) is the island ra- 
dius, F is the incoming flux, and D = i>oexTp(— Ed/ ksT) 
is the diffusion coefficient with vq being the attempt fre- 
quency, Ed the surface diffusion barrier, T the substrate 
temperature, and ks Boltzmann's constant. The solu- 
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tion is 



r 2 F 

p(r) = A + 5 In r — 



where A and B are weakly time-dependent coefficients. 
Since p(L)&p c immediately following nucleation, evalu- 
ation of A and B leads to 
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The free-adatom density on top of the islands satisfies 
Eq. (1) with a boundary condition D dp/dr | r= .R(t) = 
Sp(R(t)) where 5 = ^0 exp(— Etot/ksT) (we assume the 
same pre-factor vq for D and S). E to t = £d + £b is 
the <o<a/ barrier to move over the step edge, and Eb is 
the additional barrier (the ES barrier) which might be 
positive, zero, or even negative depending on the values 
of £ tot and E D . 7 ' 10 

The first-layer islands grow until they reach a crit- 
ical radius, R c , when second-layer free-adatom density 
reaches a critical value and nucleation on top of the first- 
layer islands occurs. Just prior to second-layer nucle- 
ation, p(0) = p c and D dp/dr \ r =R c = Sp(R c ). In this 
case the logarithmic term does not appear in the solu- 
tion of Eq. (1) and we have 
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so that from Eqs. (3) and (4) 
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From this equation, the ES barrier Eb can be 
determined. 11 

The above results are based on the assumption that 
immediately following first-layer nucleation, the inter- 
island free adatom density is given by its steady-state 
value, Eq. (2). This is unlikely and we now consider 
the consequences of relaxing this assumption. We ex- 
pect that immediately following the nucleation, p{r)ttp c 
for ro < r < L except near r = ro where p(r) is re- 
duced. However, for a large D / F ratio typical for avail- 
able experiments, 10,12-14 the adatom density quickly re- 
organizes itself with only modest island growth due to 
attachment of adatoms from the region adjacent to the 
island. We thus replace ro = 1 in Eq. (5) by a value of 
R(t) at the moment when the steady state is achieved. 
An estimate ro = l + £/10, certainly conservative given 
the dominating role of diffusion, is used below. 

We may also consider an alternative estimate of p c to 
use in place of Eq. (3). If we use the order-of-magnitude 
expression 15 p c K>L 2 F/D we obtain 
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Note that 2L can be used instead of £ in Eq. (6) (so that 
D/S = 8L 2 /R c -R c /2), cf. discussion in Ref. 15. The nu- 
(2) merical estimates based on the Eqs. (5) and (6) are given 
in Tab. I. It is worthwhile to note that the critical den- 
sity p c can in principle be determined experimentally by 
measuring the amount of material deposited before the 
onset of first-layer nucleation. It would be very useful to 
complement the observation of the onset of second-layer 
nucleation by such a measurement. Using Eq. (4), the ES 
barrier can be then estimated without any uncertainty 
connected with ro. 

To test our formulae, we employed simulations of a 
solid-on-solid KMC model in which no bulk vacancies 
or overhangs are allowed. Two processes are included 
in the model, the deposition of surface atoms and their 
migration with rates depending on their local environ- 
ment (number of lateral nearest neighbors). In the ba- 
sic version of the model, the hopping rates are given 
by an Arrhenius expression, k(E, T) = vq exp(— E/ksT), 
where E is the activation barrier to hopping comprised 
of a substrate contribution Es (always present) and a 
lateral nearest-neighbor contribution n£jv proportional 
to n, the number of the same-layer nearest neighbors in 
the initial position. 16 This model has been modified to 
include the ES barrier by considering out-of-plane next- 
nearest neighbors of an adatom before and after a hop, 17 
but for the purpose of the present investigations we used 
(4) a simpler variant in which there is an additional contri- 
bution, Eb, to the hopping barrier E for all inter-layer 
hops. The simulations were performed on 400x400 lat- 
tices with parameter values £"5 = 0.75 eV, £jv = 0.18 eV, 
i; B = 0.15 eV, at a flux £ = 0.025 monolayer(ML)/s, and 
the temperatures £ = 400 K and £ = 450 K. 

When dealing with results of simulations (and also 
with STM images) and comparing them to model pre- 
dictions, one important issue must be addressed: what 
are the proper quantities corresponding to the radius R 
and the inter-island distance 2£ of model islands for a 
collection of islands of non-ideal shapes, with a distri- 
bution of sizes and inter-island distances? After test- 
ing different possibilities, we estimated the inter-island 
distance 2£ from the island density and the critical ra- 
dius R c from the average island size. A double-check 
on the consistency of the values obtained is provided by 
the requirement that the resulting critical coverage of 
the first layer in the model given by Q C = (R C / £) 2 (note 
that an alternative C = 7ri? 2 /4£ 2 does not make an 
appreciable difference in the obtained values of the ES 
barrier) corresponds to the actual first-layer coverage at 
the onset of second-layer nucleation in a simulation or 
an experiment. 18 

The model parameters (0 C , R c , £) and the values of 
the ES barrier obtained from our simulations are sum- 
marized in Tab. I. The values of Eb extracted from the 
snapshots of the surface morphology 18 are seen to under- 
estimate (for £ = 400 K) or overestimate (for £ = 450 K) 
the value Eb = 0.15 eV actually used in simulations. 
Generally, the value of the ES barrier estimated from 
(6) Eq. (5) increases with C . At a low C , small islands are 
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TABLE I. Estimates of model parameters and the 
Ehrlich-Schwoebel barrier (in eV) using different approxima- 
tions for simulational and experimental data. The first four 
lines are results of simulations. 

C R c L Eb^ Eb° Eb c 

400 K, R, = 0.20 3.2 7.2 0.11-0.13 0.12-0.17 0.06-0.12 

400 K, R t = l d 0.30 4.3 7.8 0.10-0.12 0.11-0.16 0.05-0.11 

450 K, R t = 0.45 12.8 19.1 0.15-0.18 0.15-0.21 0.08-0.15 

450 K, R t = 2 A 0.50 14.3 20.2 0.15-0.18 0.15-0.21 0.08-0.15 

Fe(001) 12 ' 13 0.7 5.1 6.1 0.05-0.07 0.06-0.10 0.00-0.06 

Ag(lll) 10 0.55 451 608 0.18-0.22 0.19-0.23 0.13-0.19 

Pt(lll) 14 0.3 53 97 0.21-0.25 0.21-0.27 0.15-0.21 

a Eq. (5), ro = l + i/10 (first value), ro = 1 (second value). 
b Eq. (6), 2L used instead of L for the upper bound. 
c Ref. 7, Eq. (8), R = L (first value), R =2L (second value). 
d E B 1 x 0.12 - 0.14 eV (400 K) or « 0.16 - 0.19 eV (450 K) 
if the effect of incorporation is taken into account in Eq. (5), 
see text. 

responsible for the error, whereas at a high C , the on- 
set of coalescence leads to an overestimate of the barrier. 
However, we believe that the agreement (within 10-20%) 
is reasonable given the model simplifications and we did 
not attempt to compensate for these effects. 

We next look at experimental results. There is a bigger 
uncertainty in this case because an image taken exactly 
at the onset of second-layer nucleation (at a critical cov- 
erage C ) is usually not available. Fortunately, the value 
of the ES barrier is not extremely sensitive to an error 
in C and, therefore, we used the values provided by the 
authors of Refs. 7, 14 and 19. The best STM data avail- 
able are for Fe/Fe(001) homoepitaxy. 12,13 We used the 
results for growth of Fe at 20° C and estimated the av- 
erage inter-island distance 2L from the observed density 
of islands k, 8xl0 12 cm -2 (Ref. 12). The critical radius 
R c was estimated from the first-layer coverage 19 of 0.7 
ML at a total coverage of 0.8 ML. The value of the ES 
barrier we obtained (Tab. I) is close to the estimates of 

0.05 eV based on comparisons between the evolution of 
the surface roughness for the same STM data and KMC 
simulations. 20 Our estimate may be too high because C 
is rather high. Other available experimental data are for 
Pt(lll) from Ref. 14 and for Ag(lll) from Refs. 7 and 
10. We used the estimates for L and C provided by 
the authors of these publications and our estimates of 
R c from the critical coverage and their STM pictures. 
The results are summarized in Tab. I. 21 The values of 
© c suggest that the ES barrier for Pt(lll) may in fact 
be even higher, and the one for Ag(lll) lower than our 
results. 

When extracting values for physical constants based 
on comparison with a simple model (analytical or KMC), 
it is always important to try to estimate how sensitive 
they are to the model simplifications. For real experi- 
mental systems, additional factors not considered in the 
model used here might play a role and it seems most 
of these lead to underestimating the ES barrier. Let us 
briefly discuss some of the possibilities: 



(i) A "leaking" ES barrier. It has been proposed based 
on microscopic calculations that the barrier to descend 
might be lower 22 or higher 23 for rough steps or near 
kinks. Also, there are two types of steps (usually called 
A and B steps) on fcc(lll) surfaces and it has been found 
in Refs. 22-24 that the barrier is much lower for B-type 
steps. 25 The value of the ES barrier obtained from a sim- 
ple model is an effective value (averaged over the edges 
of all islands) and may be thus rather different from a 
barrier at a straight (e.g., A-type) step. 

(ii) Fractal islands. If the island morphology is signifi- 
cantly different from the compact one we assumed, the 
model cannot be applied straightforwardly. Meyer et al. 7 
suggested that in the case of fractal morphology the is- 
land radius can be replaced with the width of fractal 
arms but this is certainly a very crude approximation. 
In any case, the barrier obtained from the model as out- 
lined above will underestimate the real barrier of fractal 
islands because of more frequent attempts by adatoms 
to hop down the steps. A recently suggested effect, 22 
that an adatom can be trapped near a descending step, 
will also lead to more frequent attempts to hop down 
the step and thus a too low value of the barrier obtained 
from the model used here. 

(iii) Post-deposition non-diffusional incorporation. In 
some cases, some of the atoms deposited near island 
edges can descend to a lower layer directly following their 
impact onto the surface. 22,26,27 This effect will again lead 
to a smaller barrier obtained from the model as com- 
pared to the real one, but the error will depend strongly 
on the average island size. We can provide a quantita- 
tive estimate by modifying the model such that all atoms 
deposited near island edges within a distance Ri (in the 
region R — Ri <r < R, see Fig. 1) immediately descend 
to the lower layer instead of staying on top of the island. 
We then have to replace the flux F on top of the island 
by F', F' = F(l - 2Ri/R + Rf/R 2 ). This does not af- 
fect Eq. (3), but in Eq. (4), F must be replaced by F' 
so that all but the last term on the right-hand side of 
Eq. (5) are multiplied by the ratio F/F'. To test the 
effect of incorporation, we repeated our simulations with 
a different deposition rule: instead of simply dropping 
a particle onto a randomly selected site, we search for a 
site with the highest coordination and at a lower layer 
within a square region of a linear dimension of 2i? 8 - + 1 
(Ri =1 and 2 were used for the temperature 400 K and 
450 K, respectively). The results given in Tab. I show 
that the incorporation leads to a change in the estimate 
of the barrier of the order of 0.02 eV (400 K) and 0.01 eV 
(450 K). Note however, that even incorporation over the 
distance of several lattice constants would not change 
appreciably the estimates of the ES barrier for Ag(lll) 
and Pt(lll) because large island sizes make this effect 
unimportant. 

In conclusion, we have investigated a simple model to 
determine the value of the activation barrier to interlayer 
transport from information on the surface morphology 
in the initial stages of homoepitaxial growth. The model 
has been tested using kinetic Monte Carlo simulations 
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of homoepitaxial growth and found to provide good esti- 
mates of the value of the step-edge barrier. This method 
of extracting the step-edge barrier from experimental 
data is straightforward and easy to use: the only in- 
formation necessary is the island density and first-layer 
(or total) coverage close to the onset of nucleation on 
top of first-layer islands. Although we used scanning- 
tunneling microscopy data, high-resolution diffraction 
measurements can be used as well. 28 Various processes 
not accounted for in the simple model we use lead usually 
to underestimating the real value of the barrier. 
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